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Abstract 

The hemoglobin of jawed vertebrates is a heterotetramer protein that contains two a- and two p-chains, which are encoded by 
members of oc- and p-globin gene families. Given the hemoglobin role in mediating an adaptive response to chronic hypoxia, it is likely 
that this molecule may have experienced a selective pressure during the evolution of cetaceans, which have to deal with hypoxia 
tolerance during prolonged diving. This selective pressure could have generated a complex history of gene turnover in these clusters 
and/or changes in protein structure themselves. Accordingly, we aimed to characterize the genomic organization of a- and p-globin 
gene clusters in two cetacean species and to detect a possible role of positive selection on them using a phylogenetic framework. 
Maximum likelihood and Bayesian phylogeny reconstructions revealed that both cetacean species had retained a similar complement 
of putatively functional genes. For the oc-globin gene cluster, the killer whale presents a complement of genes composed of HBZ, HBK, 
and two functional copies of HBA and HBQ genes, whereas the dolphin possesses HBZ, HBK, HBA and HBQ genes, and one HBA 
pseudogene. For the p-globin gene cluster, both species retained a complement of four genes, two early expressed genes — HBE and 
HBH — and two adult expressed genes — HBD and HBB. Our natural selection analysis detected two positively selected sites in the HBB 
gene (56 and 62) and four in HBA (1 5, 2 1 , 49, 1 20). Interestingly, only the genes that are expressed during the adulthood showed the 
signature of positive selection. 
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Introduction 

Hemoglobin (Hb) is one of the best studied of all macromo- 
lecular proteins. The main function of Hb is to transport mo- 
lecular 0 2 from the respiratory systems to the tissues, where it 
is released for cell use (Dickerson and Geis 1983). The Hb of 
jawed vertebrates is a heterotetramer protein that contains 
two oc- and two p-chains, which are encoded by members 
of the corresponding a- and p-globin gene families. These 
families are ontogenetically regulated, which means that func- 
tionally distinct Hb isoforms are differentially expressed in em- 
bryonic and adult erythroid cells. They are also biochemically 
optimized for oxygen transport under different physiological 
conditions faced during all stages of development (Hardison 
1998; Storz et al. 2008; Hoffmann et al. 2010). During the 
course of vertebrate evolution, the oc- and p-globin gene 



families have undergone a complex history of gene duplica- 
tion and divergence that have given rise to a great diversity in 
extant species (Hoffmann et al. 2008, 2010; Opazo et al. 
2008a,b). These regulatory and functional divergences consti- 
tute an important source of variation that affects 
physiologically important aspects of hemoglobin properties 
(Hoffmann et al. 2008; Storz et al. 2008). From a genomic 
standpoint, oc- and p-globin gene families are located in two 
different clusters in the genomes of amniotes. 

Physiological adaptations to hypoxia are necessary in en- 
vironments where oxygen availability is reduced. For exam- 
ple, evidence from a number of high-altitude vertebrates 
indicates that modification of Hb properties typically plays 
a key role in mediating an adaptive response to chronic hyp- 
oxia: high-altitude species generally have Hbs with higher 0 2 
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binding affinities than those of their lowland relatives 
(e.g., Weber et al. 2002; Storz 2007). In aquatic environ- 
ment, diving vertebrates also have developed suitable mech- 
anisms for the maintenance of an adequate 0 2 supply to 
tissues under hypoxic conditions (Tellone et al. 2000). For 
example, Hbs from emperor penguin (Meir and Ponganis 
2009) and sea turtle (Petruzzelli et al. 1996) display a Bohr 
effect that appears well adapted for gas exchange during 
prolonged dives. It is likely that the Hb of diving marine mam- 
mals, which usually show a remarkable hypoxia tolerance, 
would possess adaptive changes as well. In this respect, 
cetaceans may have developed particular molecular mecha- 
nisms for the maintenance of adequate oxygen supply to 
tissues during acute hypoxia (Brix et al. 1990). 

Although a few cetacean a- and p-globins have been 
sequenced so far and analyzed functionally (Brix et al. 1990; 
Tellone et al. 2000; Corda et al. 2003; Remington et al. 2007; 
Manconi et al. 2009), no previous work has characterized the 
a- and p-globins clusters on this particular lineage. Most of the 
earlier studies have shown that the oxygen-binding behavior of 
cetaceans' hemoglobin has unusual properties probably linked 
to their peculiar diving behavior: lower intrinsic oxygen affini- 
ties and greater Bohr effects. These characteristics could en- 
hance oxygen release to respiring tissues when exercise causes 
a drop in levels of ATP. Remington etal. (2007) highlighted that 
the relative low-oxygen affinities are intrinsic to the Hb's struc- 
ture and showed that genetic diversity was clearly responsible 
for the variations in Hb type and function of these animals. 

In this context, it is tempting to think that Hb may have 
experienced a greater selection pressure to modify its func- 
tional properties, resulting in a molecule more adapted to 
prolonged diving behavior. Thus, oc- and p-globin gene fami- 
lies are great candidates to determine whether the unusual 
properties of cetacean's hemoglobin are attributable to differ- 
ences in Hb gene repertoire and/or to positive selection in 
amino acid residues that play important roles in controlling 
Hb-0 2 affinity (e.g., Perutz 1983; Berenbrink 2006; Storz 
and Moriyama 2008; Mandic et al. 2009). Accordingly, the 
objectives of this study were as follows: 1) to characterize the 
organization of the Hb gene clusters in cetaceans, 2) to infer 
orthologous relationships among duplicated copies of a- and 
P-globin genes, and 3) to infer the potential role of positive 
selection on the a- and p-globin genes in the cetacean lineage 
using a phylogenetic framework. 

Materials and Methods 

Sequence Data 

We used bioinformatic tools to identify the full complement of 
structural genes in the oc- and p-globin gene families in 13 
laurasiatherian species and one primate as outgroup species 
(supplementary table S1, Supplementary Material online). We 
included four cetartiodactyls (dolphin, Tursiops truncatus; 



killer whale, Orcinus orca; cow, Bos Taurus; pig, Sus scrota); 
two bats (megabat, Pteropus vampyrus; microbat, Myotis luci- 
fugus); four carnivores (cat, Felis catus; dog, Canis familiaris; 
panda, Ailuropoda melanoleuca; ferret, Mustela putorius); 
one perissodactyl (horse, Equus caballus); and two insectivores 
(hedgehog, Erinaceus europaeus; shrew Sorex araneus). 
Genes were manually annotated by comparing known exon 
sequences to genomic contigs using the program Blast2seq, 
version 2.2 (Tatusova and Madden 1999), and also using the 
program Genscan (Burge and Karlin 1997). All sequences 
were aligned using the program MUSCLE (Multiple 
Sequence Comparison by Log-Expectation) (Edgar 2004). 
Putatively functional genes were characterized by an intact 
open reading frame with the canonical three exon/two 
intron structure typical of vertebrate globin genes, whereas 
pseudogenes were identifiable because of their high sequence 
similarity to functional orthologs and the presence of inacti- 
vating mutations and/or the lack of exons. 

Inferring Orthologous Relationships 

We inferred orthologous relationships by reconstructing 
phylogenetic relationships among the a- and p-globin genes 
using Bayesian and maximum likelihood approaches, as im- 
plemented in Mr.Bayes v.3.1.2 (Ronquist and Huelsenbeck 

2003) and Treefinder version October 2008 (Jobb et al. 

2004) . The models of nucleotide substitution were selected 
through the "propose model" tool of Treefinder version 
October 2008 (Jobb et al. 2004). For the Bayesian analyses, 
two simultaneous independent runs were performed for 
30,000,000 iterations of a Markov chain Monte Carlo algo- 
rithm, with six simultaneous chains, sampling every 1,000 
generations. Support for the nodes and parameter estimates 
were derived from a majority rule consensus of the last 1 5,000 
trees sampled after convergence. In maximum likelihood, we 
estimated the best tree using the models of nucleotide sub- 
stitution previously selected, and support for the nodes was 
estimated with 1,000 bootstrap pseudoreplicates. Because of 
the occurrence of interparalog gene conversion among globin 
genes (Storz et al. 2007; Opazo et al. 2008b), we used phy- 
logeny reconstructions of noncoding sequences (intron 2 and 
flanking sequences, 1 kb upstream of the start codon, and 
1 kb downstream immediately after the stop codon) to infer 
orthologous relationships among oc- and p-globin genes. 

Detecting Positive Selection 

To investigate the possible role of natural selection in the evo- 
lution of oc- and p-globin genes in the cetacean lineage, we 
used the branch-site model as implemented in the program 
codeml included in the software package PAML v4.4 (Yang 
2007). Because the branch-site analysis estimates rates of 
evolution on a codon-by-codon basis on a specific branch of 
the tree, its implementation is particularly useful in cases when 
different gene segments evolve at different rates. This model 
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Fig. 1. — Genomic structure of a- and p-globin gene clusters in laurasiatherian species. The orientation of the cluster is from 5' (on the left) to 3' (on the 
right). Diagonal slashed on the beta-globin cluster of the killer whale indicates gaps in genomic coverage. Because of these gaps, the exon 3 of the HBD gene 
of the killer whale is missing. 



assumes that branches on the tree are divided into foreground 
branches, where some sites may be under positive selection, 
and background branches where positive selection is absent 
(Yang and Nielsen 2002; Yang et al. 2005; Zhang et al. 2005). 
Under this methodology, sites are categorized into four clas- 
ses 0, 1 , 2a, and 2b with proportions of pO, p1 , p2a, and p2b, 
respectively. Site class 0 includes codons that evolve under 
purifying selection on both the foreground and background 
branches, with 0 < co 0 < 1 . In site class 1 , codons evolve neu- 
trally in all lineages, with coi = 1 , whereas in classes 2a and 2b, 
positive selection is allowed on the foreground branches with 
co 2 > 1, but not on the background branches. This model is 
compared with the corresponding null hypothesis of neutral 
evolution, where co 2 is fixed to 1 . If the null hypothesis is re- 
jected by the likelihood ratio test (LRT), a Bayes empirical Bayes 
approach is used to calculate the posterior probabilities that 
each site has evolved under positive selection on the fore- 
ground lineage (Yang et al. 2005). In all cases, three starting 
co values (0.5, 1, and 2) were used to check the existence of 
multiple local optima. In our case, the ancestral branch of the 
cetacean clade was labeled as the foreground branch. All 



the analyses were based on a phylogeny that includes repre- 
sentative species of laurasiatherian mammals. The tree topol- 
ogy used to conduct the analyses of variable co among 
lineages and sites is based on published literature (Nery 
et al. 2012). 

Results 

Genomic Structure of the Dolphin and Killer Whale 
a- and (3-Globin Gene Clusters 

The oc-globin gene cluster of the cetacean species included in 
our study appeared to have retained a similar complement of 
putatively functional genes: the killer whale presents a com- 
plement of genes composed of HBZ, HBK, and two func- 
tional copies of HBA and HBQ globin genes, whereas the 
dolphin possesses HBZ, HBK, HBA and HBQ genes, and 
one HBA pseudogene (fig. 1). Phylogenies based on noncod- 
ing sequences indicate that the oc-globin genes found in the 
dolphin and in the killer whale are 1:1 orthologs to the oc- 
globin gene repertoire present in other mammals (fig. 2). The 
oc-globin gene cluster of the cetaceans retained very similar 
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Fig. 2. — Maximum likelihood trees showing relationships among a-globin genes based on 1 kb of 5 ; flanking sequences (left), intron 2 (center), and 1 kb 
of 3' flanking sequence (right). Values on relevant nodes denote bootstrap support values and Bayesian posterior probabilities. 



2362 Genome Biol. Evol. 5(1 2):2359-2367. doi:1 0.1 093/gbe/evt1 76 Advance Access publication November 20, 201 3 



Globin Evolution in Two Cetacean Lineages 



GBE 



patterns of intergenic spacing in comparison to other laura- 
siatherian species: from the initiation codon of HBZ to the 
termination codon of HBQ, the gene cluster covered 
11 ,690 bp in the dolphin, 10,601 in the killer whale, and 
12, 159 bp in the pig. Among laurasiatherian mammals, the 
number of putatively functional genes ranged from three 
in the ferret to eight in the cat. All species included in this 
study possess at least one functional copy of HBZ and HBA 
(fig. D. 

For the p-globin gene cluster, cetaceans retained a com- 
plement of four putatively functional genes (from 5' to 30, 
two early expressed genes — HBE and HBH — and two adult- 
expressed genes — HBD and HBB (fig. 1). Phylogenies based 
on noncoding sequences indicate that the p-globin genes 
found in both the dolphin and the killer whale are 1:1 ortho- 
logs to the p-globin gene repertoire in other mammals (fig. 
3). Also the dolphin p-globin gene cluster has remained sim- 
ilar in terms of intergenic spacing in comparison to other 
phylogenetically related species. From the first codon of 
HBE to the termination codon of HBB, the dolphin cluster 
spans 20,317 bp and the pig 19,51 6 bp. Unfortunately, the 
current state of the killer whale genome assembly does not 
permit inferences regarding patterns of intergenic spacing in 
the p-globin gene cluster. Among laurasiatherians, gene 
repertoire of the p-globin gene cluster varies from two 
functional genes in the pig up to eight in the cow, whose 
globins were involved in en bloc duplications (Townes et al. 
1984). 

Variation in Omega Ratio 

The results of the branch-site analyses are summarized in 
table 1. It is interesting to highlight that for both adult- 
expressed genes (HBA and HBB) the LRT indicates that 
the models which estimate a class of sites with an co value 
higher than 1 had a significantly better fit than the null 
model in which the co value was fixed to 1 
[2A€ = 2*(^ -l 0 )= 11.11 and 9.49; p = 0.008 and 0.002, 
respectively]. In these analyses, we were able to identify sites 
under positive selection in both genes. In the case of the HBA 
gene, sites 15, 21, 49, and 120 were positively selected, 
whereas for the HBB, gene sites 56 and 62 in HBB were in- 
ferred under the action of positive selection. In the case of the 
embryonically expressed genes (HBZ and HBE), the model A 
did not fit the data better than the null model in which the co 
value was fixed to 1 [2 Al = 2*(l< [ -l 0 ) = 0 and 0.08; p = 1 .0 
and 0.69, respectively]. 

Discussion 

The transition from land to an aquatic environment in the early 
evolutionary history of cetaceans involved extensive modifica- 
tions on their physiological, anatomical, and behavioral sys- 
tems (Butler 2004; Ramirez et al. 2007). One of the most 
challenging aspects of a breath-holding animal living in an 



aquatic environment is to deal with extended periods of lim- 
ited oxygen during submersion. Cetaceans have developed 
several strategies to cope with this limitation, such as 0 2 
saving by reduction of the metabolic rate, selective vasocon- 
striction to assure 0 2 supply in sensitive organs, massive 
overexpression of myoglobin in aerobic muscles, a modified 
hemoglobin that performs better under acute hypoxia condi- 
tions, among others (Brix et al. 1990; Ramirez et al. 2007). 
Given that modifications of Hb function are often implicated 
in adaptation to acute hypoxia and because much is known 
about the structure-function relationships of vertebrates Hbs 
and their role in blood-0 2 transport, the study of Hb function 
in vertebrate species that have developed great hypoxia toler- 
ance, such as cetaceans, provides a great opportunity to elu- 
cidate detailed molecular mechanisms of physiological 
adaptation. 

Cluster Organization in the Cetacean Lineage 

The laurasiatherian globin gene repertoire typically comprises 
four functional genes for the oc- (HBZ, HBK, HBA, HBQ) and 
P-globin gene clusters (HBE, HBG, HBB, HBD). In this regard, 
the cetacean species included in our study, in contrast to other 
lineages that underwent a complex history of gene turnover, 
present a typical gene repertoire for a laurasiatherian species 
(fig. 1). Not surprisingly, for the dolphin and the killer whale 
oc- and p-globin clusters, the embryonic globin genes were 
located at 5 r of those with an adult expression, as is the gen- 
eral pattern in the globin gene clusters of vertebrates. 

It has been suggested that the variation in the globin gene 
repertoire among different lineages may constitute an impor- 
tant source of variation that affects physiologically important 
aspects of blood oxygen transport (Weber 2007; Storz and 
Moriyama 2008; Storz et al. 2008). As diving mammals, which 
have to face prolonged periods of hypoxia during diving, one 
could expect that the globin genes in the cetacean lineage 
would show a complex history of gene turnover. But instead 
they appear to have retained a typical gene repertoire as other 
laurasiatherian mammals (fig. 1). Given this pattern, it is likely 
that the phenotypic differences are most attributable to 
changes in protein structure rather than in copy number 
variation and genomic organization. 

Variation in Selective Pressure 

The oc- and p-globin gene clusters are ontogenetically regu- 
lated and biochemically optimized for oxygen transport under 
the different physiological conditions that are encountered 
during the embryonic and adult stages of development. It is 
well known that the fetal hemoglobin is structurally different 
from the adult hemoglobin by possessing greater affinity for 
oxygen than the adult hemoglobin (Delivoria-Papadopoulos 
and McGowan 1998). As a consequence, fetal hemoglobin 
combines more rapidly with oxygen at low tension than does 
adult hemoglobin. This is what is needed because the partial 
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Fig. 3. — Maximum likelihood trees showing relationships among p-globin genes based on 1 kb of 5' flanking sequences (left), intron 2 (center), and 1 kb 
of 3' flanking sequence (right). Values on relevant nodes denote bootstrap support values and Bayesian posterior probabilities. The absence of killer whale 
HBB on the "1 kb upstream" tree is due to the poor quality of its genomic sequence. 



pressure of oxygen in the arterial blood is considerably lower 
than that of the atmospheric environment at the sea level. It is 
then expected that both the embryonic and the adult globins 
might be subject to different selective regimes. Our selective 
pressure analyses reflected exactly this scenario: only the ce- 
tacean adult a- (HBA) and p-globin (HBB) presented signatures 



of positive selection, whereas the embryonic a- (HBZ) and p- 
globins (HBE) seem to be more constrained (table 1). This 
evolutionary pattern (i.e., embryonic genes evolving more 
slowly than those expressed in adults) is common in the liter- 
ature, and it seems to be widespread on vertebrate genome 
(Goodman 1963; Shapiro 1991; Roux and Robinson-Rechavi 
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Table 1 

Log Likelihood Values and Parameters Estimates under Different Branch-Site Models, Where the Cetacean Lineage Was Labeled as Foreground 
Branch 



Gene 



In L 



Model 



Estimated parameters 



p 0 = 0.760, Pt = 0.238 p 2a = 0, 

p 2b = 0.001; coo = 0.001, co 2a = 1-000, co 2b = 1-000 
p 0 = 0.761 , = 0.239, p 2a = 0, p 2b = 0; 

co 0 = 0.064, co 2a = 1 .000, co 2b = 1 .000 



Positively Selected 
Sites 



HBZ (embryonic) 



-2490.855 
-2490.855 



Null 

Model A 



p 0 = 0.778, Pt = 0.191, p 2a = 0.024, p 2b = 0.006; 

coo = 0.050, cot = 1.000, co 2a = 1.000, co 2b = 1.000 
p 0 = 0.778, p 1 = 0.184, p 2a = 0.031, p 2b = 0.007; 

(Oo = 0.050, 03! = 1.000, 03 2a = 28.565, o» 2b = 28.565 



NA 
NA 



HBA (adult) 



-2327.998 
-2322.444 



Null 

Model A 



p 0 = 0.797, Pt = 0.166, p 2a = 0.031, p 2b = 0.006; 

co 0 = 0.049, co-| = 1 .000, co 2a = 1 .000, co 2b = 1 .000 
p 0 = 0.823, Pt = 0.171, p 2a = 0.005, p 2b = 0.001; 

o = 0.049, a.! = 1 .000, 03 2a = 1 9.443, 03 2b = 1 9.443 



NA 

15, 21, 49, 120 



HBE (embryonic) 



-2003.994 
-2003.917 



Null 

Model A 



03 0 = U.U4y, 03t = I.UUU, 03 2a = \y.*V\S, 03 2b = iy/w. 

p 0 = 0.651, Pt = 0.263, p 2a = 0.060, p 2b = 0.024; 

03 0 = 0.053, 03! = 1 .000, 03 2a = 1 .000, 03 2b = 1 .000 
p 0 = 0.702, p 1 = 0.278, p 2a = 0.014, p 2b = 0.005; 

O3o = 0.057, 031 = 1.000, 03 2a = 24.033, 03 2b = 24.033 



NA 
17 



HBB (adult) 



-2296.787 
-2292.040 



Null 

Model A 



NA 

56, 62 



note. — In L likelihood value, p proportion of site class, co omega value for site classes. Sites inferred under positive selection for HBA gene had posterior probabilities 
values of 0.60, 0.88, 0.56, and 0.55, respectively. For HBB gene, posterior probabilities values were 0.79 and 0.99, respectively. 



2008). There are two nonmutually exclusive hypotheses to 
explain this pattern. The "developmental constraint" hypoth- 
esis suggests that genes expressed in early developmental pro- 
cesses are under strong negative selection as a means to avoid 
mutations that could cause deleterious cascading effects (Raff 
1996). The "selection opportunity" hypothesis suggests that 
later stages in ontogeny give rise to greater opportunity for 
selection to act due to the exposure to varying environments 
(Darwin 1872; Gould 1977). Our results support both predic- 
tions of these models (embryonic globins under negative 
selection and adult globins subject to greater positive selection 
pressure). 

The branch-site analyses identified four positively selected 
sites in the cetacean HBA gene— 15 (Ser), 21 (Ser), 49 (Gly), 
and 120 (Ser) — and two sites in the cetacean HBB gene — 56 
(Lys) and 62 (Lys). Interestingly, almost all these substitutions 
are shared by both cetacean species, indicating that positive 
selection in the common ancestor of these species was 
responsible for the remodeling of this protein (fig. 4). It is 
known that Hb-oxygen affinity can be modified by amino 
acid substitutions that decrease or enhance the structure for 
oxygen binding or by changing the affinity of Hb for allosteric 
effectors (Perutz 1989). The residue 62 is a positively selected 
site in HBB gene and is known to be part of the E helix of the 
beta chain. According to the three-dimensional model of 
Perutz et al. (1968), some residues of the E helix seem to be 
essential in the maintenance of the heme in the nonpolar 
pockets of the alpha and beta chains. 

Regarding the positively selected sites in the HBA gene, the 
residue 21 (Ser) is almost always occupied by an alanine or a 



valine on eukaryotes, and in few species by leucine or isoleu- 
cine (Harteveld et al. 2007). The presence of a serine at this 
position found in both cetacean species (fig. 4) was already 
described as naturally occurring in human but the physiolog- 
ical consequences for this specific mutation are still unknown 
(Harteveld et al. 2007). The site 120 (Ser) is known to be in- 
volved in the ocipi contacts in the Hb molecule (Perutz and 
Lehmann 1968; Sacks et al. 1978), and consequently, an 
amino acid substitution at this site causes abnormality for 
oxygenation, as already described by studies that found new 
variants at this residue (Harano et al. 1 989; Dingol et al. 2006). 

We did not obtain the oxygen dissociation curve from the 
dolphin Hb and the physiological effects of the positively 
selected residues 56 in HBB, and 15 and 49 in HBA are yet 
to be determined. However, other studies already described 
natural mutants occurring at these positions in human, with 
different amino acids (Szelenyi et al. 1980; Turbpaiboon et al. 
2002; Chang et al. 2002; Akar et al. 2003; Williams et al. 
2007). 

Although not identified as a positively selected site, we 
noted that the residue 65 of the HBB gene is occupied by a 
glutamine residue (Gin) only in the dolphin species (fig. 4). This 
residue is located very close to the distal histidine E7, an in- 
variant amino acid site, which is implicated in the oxygenation 
process (Rohlfs et al. 1990). Garel et al. (1976) described a 
human hemoglobin variant in which a glutamine residue — the 
same found in the dolphin — substitutes the lysine residue in 
position 65. They found that this new variant results in a mod- 
erate decrease in cooperativity without changing the Hb sta- 
bility. It suggests that hemoglobin with a 65Gln will result in a 
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t t 2 i° 3 i° 4 .° 5 .° 6 .° 7 .° 8 .° 9 i° 1 ?° 4° ^ f i f 

Cow HBZ SLTRAERTIIVSMWSKISTQADLIGTETLERLFSCYPQAKTYFPHFDLHTGSAQLRAHGSRVVAAVGDAVKSIDNVTSALSKLSELHAYVLRVDPVNFKFLSHCLLVTLASHFPADFTADVHAAWDKFLSIVSGVLTEKYR 

Dolphin HBZ G I A S...T A K AG L V A R 

Killer whale HBZ V . SP . DK . NVKGT . A. . GNHSAEY . A. A. . .M. INF. ST S GH IKG. . KK.AD.LTK. .GH. . . LPD . . . E . .D. . . HK L L.L PS . . . SL AS. .T. . .S. . . 

Cow HBA VLSAADKGNVKAAWGKVGGHAAEYGAEALERMFLSFPTTKTYFPHFDLSHGSAQVKGHGAKVAAALTKAVEHLDDLPGALSELSDLHAHKLRVDPVNFKLLSHSLLVTLASHLPSDFTPAVHASLDKFLANVSTVLTSKYR 

Dolphin HBA ...P...T...GT.S.I.N.S IN..S S G I K...D G.I.N..D C L...A S S 

Killer whale HBA ...P...T...GT.A.I.N.S IN..S S G I K...D G.I.N..D C L...A S S 

1 if 20 30 40 50 60 70 80 90 100 1J0 IfO 130 140 

Cow HBE VHFTGEEKAAITDLWGKMNMEEAGGEALGRLLVFYPWTQRFFDSFGNLSSASAIMGNPKVKAHGKKVLTSFGEAMKNSDNLKGAFAKLSELHCDKLHGDPENFRLLGDVIVIVLATHSGREFTPDMQAAWQKLVSGVATALAHKYH 

Dolphin HBE. ...A G...RV.V V N P D.V..M V N.L..I..S.F EL A 

Killer whale HBE A G...RV.V V N P D.V..M V N.M..I..S.F EL A 

Cow HBB VHFTAEEKAAITGLWGKVNVEEAGGEALGRLLVVYPWTQRFFDSFGNLSSASAIMGNPKVKAHGKKVLTSFGEAIKNLDNLKGAFAKLSELHCDKLHVDPENFRLLGNVIVIILATHFGREFTPDVQAAWQKLVSGVATALAHKYH 

Dolphin HBB..L.G...S.V.A V E...D..T.D.V.K..N..IK1. .Q...A GL.H..D...T..A L.VV..R...K EL.S.Y..V.A 

Killer whale HBB..L.G...S.V.A V E...D..T.D.V.K..N..|k| A GL.H..D...T..A L.VV..R...K EL.S.Y..V.A 



Fig. 4. — An alignment of HBA and HBB amino acid sequences from cow, dolphin, and killer whale. Positively selected amino acid sites are shaded. 



protein with decreased oxygen affinity. In diving birds and 
mammals, the Hb-0 2 affinities are not high, because they 
do not experience environmental hypoxia while breathing 
(i.e., low 0 2 partial pressure; Snyder 1983; Willford et al. 
1 990; Ramirez et al. 2007). These animals may instead benefit 
from an Hb with lower intrinsic oxygen affinities, which would 
enhance Hb-0 2 unloading during acute hypoxia, and also a 
high Bohr effect, which would intensify 0 2 unloading when 
acidosis increases during diving (Willford et al. 1990; Butler 
and Jones 1997; Ramirez et al. 2007). 

It is important to note that inferences about functional 
divergence based on omega variation among lineages do 
not replace research on the biochemical properties of proteins. 
Nevertheless, our analyses is the first informative step toward 
understanding the functional divergence and produce candi- 
dates to perform site-directed mutagenesis studies and con- 
firm whether or not these specific amino acid substitutions 
confer cetacean's hemoglobin with unique properties that 
underlie their adaptation to the acute hypoxia. 

General Conclusions 

Taken together, our analyses represent a step toward increas- 
ing our knowledge about the oc- and p-globin genes in the 
cetacean lineage and how they evolved. Our study revealed 
that oc- and p-globin gene families of the dolphin and killer 
whale do not present a complex history of duplication and 
divergence; instead, they appeared to have retained a typical 
laurasiatherian gene repertoire. Our branch-site analyses were 
able to identify positively selected sites only across the adult 
oc- and p-globin genes, reinforcing the interesting pattern 
already described for other species, where embryonic genes 
appear to evolve slower than the adult expression genes 
(Goodman 1963). 

Supplementary Material 

Supplementary table S1 is available at Genome Biology and 
Evolution online (http://www.gbe.oxfordjournals.org/). 
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